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We propose to utilize multiple-graphene-layer structures with lateral p-i-n junctions for terahertz 
(THz) and infrared (IR) photodetection and substantiate the operation of photodetectors based on 
these structures. Using the developed device model, we calculate the detector dc responsivity and 
detectivity as functions of the number of graphene layers and geometrical parameters and show that 
the dc responsivity and detectivity can be fairly large, particularly, at the lower end of the THz 
range at room temperatures. Due to relatively high quantum efficiency and low thermogeneration 
rate, the photodetectors under consideration can substantially surpass other THz and IR detectors. 
Calculations of the detector responsivity as a function of modulation frequency of THz and IR radia- 
tion demonstrate that the proposed photodetectors are very fast and can operate at the modulation 
frequency of several tens of GHz. 
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I. INTRODUCTION 

Unique properties of graphene layers (GLs) P, [H, Q 
make them promising for different nanoelectronic device 
applications. The gapless energy spectrum of GLs, which 
is an obstacle for creating transistor-based digital cir- 
cuits, opens up prospects to use GLs in terahertz (THz) 
and infrared (IR) devices. Novel optoelectronic THz and 
IR devices were pro pos ed and evaluated, in particular, in 
Refs. 0, la fi 0, la H, LLIJ, UJJ ■ Recent success in fabricating 
multiple-GL structures with long momentum relaxation 
time of electrons and holes [12j | promises a significant 
enhancement of the per formance of futures graphene op- 
toelectronic devices [l3j . 

In this paper, we study the operation of THz and IR 
photodetectors based on multiple-GL structures with p- 
i-n junctions. We refer to the photodetectors in ques- 
tion as to GL-photodetectors (GLPDs). We focus on 
GLPDs with p- and n-doped sections in GLs near the 
side contacts [14], HH , p + and n + contacts (for example, 
made of doped poly-Si [h|). and multiple GL-structures 
with the Ohmic side contacts and split-gates which pro- 
vide the formation of the electrically induced p- and n- 
sections [H GJ, El US HH [Hj] ■ The device structures 
under considerations are shown in Figs. 1(a) and 1(b). 
It is assumed that the highly-conducting GL(s) between 
the SiC substrate and the top GLs is removed. Mul- 
tiple GL-structures without this highly conducting GL 
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can be fabricated using chemical/mechanical reactions 
and transferred substrate techniques (chemically etching 
the substrate and the highly conducting bottom GL [2j| 
or mechanically peeling the upper GLs, then transferring 
the upper portion of the multiple-GL structure on a Si 
substrate) . 

Using the developed device model, we calculate the 
GLPD responsivity and detectivity as THz or IR pho- 
todetector (Sees. II and III), evaluate its dynamic re- 
sponse (Sec. IV), and compare GLPDs with some other 
THz and IR photodetectors (Sec. V), in particular, 
with quantum-well infrared photodetectors (QWIPs) 
and quantum-dot infrared photodetectors (QDIPs). In 
Sec. VI, we discuss possible role of the Pauli blocking 
in the spectral characteristics of GLPDs and draw main 
conclusions. The calculations related to the effect of 
screening of the vertical electric field on the formation 
of the electrically induced p- and n-sections in multiple- 
GL structures are singled out in the Appendix. 



II. RESPONSIVITY AND DETECTIVITY 



We assume that the intensity of the incident THz 
or IR radiation with the frequency Q apart from the 
dc component Iq includes the ac component: I(t) = 
Iq + SI U exp(— iu>)t, where 51^ and w are the amplitude 
of the latter component and its modulation frequency, 
respectively. In such a situation, the net dc current (per 
unit width of the device in the direction perpendicular to 
the current) can be presented in the following form: 
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FIG. 1: Device structures of GLPDs with (a) doped and (b) 
electrically induced p — i — n- junctions. 
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If if = 1, Eq. (5) at fl/2iT = 1 THz, yields, R ~ 
12 A/W. For K = 50 - 100 (if* ~ 30 - 39), setting 
if = 50, in the frequency range S7/27T =1 — 10 THz, 
from Eq.(4) we obtain R ~ 35 - 350 A/W. 

The dark-current limited detectivity, Z)*, defined as 
£>* = (JP hoto /NP)y/A ■ A/, where A is the noise (in 
Amperes), P is the power received by the photodetec- 
tor (in Watts), A is the area of the photodetector (in 
cm 2 )(see, for instance, [H, Hi|), and Af is the band- 
width, can be expressed via the responsivity i? and the 
dark current jQ ark H, where H is the device width in the 
direction perpendicular to the current, can be presented 
as 
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Here e is the electron charge, K is the number of GLs, 
21 is the length of the GL's i-section, g t h and g tU nn = 
(eV/2lv^ 3 fi) 3 / 2 /87r 2 are the rates of thermal and tun- 
neling generation of the electron- hole pairs (per unit 
area) [22|, V is the bias voltage, vw — 10 8 cm/s is the 
characteristic velocity of electrons and holes in graphene, 
fin = @W — 2f(Ml/2)] is the absorption coefficient of ra- 
diation in a GL due to the intcrband transitions [24| . 
where fi — (ire 2 /ch) ~ 0.023, f(e) is the distribution 
function of electrons and holes in the i-section, H is the 
Planck constant, and c is the speed of light. The term 
gtunn in Eq. (2) is associated with the interband tunnel- 
ing in the electric field in the i-section. The quantity 

= J (l - fi) K - k is the intensity of THz or IR radia- 
tion at the fc-th GL (1 < k < K). 

At a sufficiently strong reverse bias, electrons and holes 
are effectively swept out from the i-section to the contacts 
and heated. Due to this, one can assume that under the 
GLPD operation conditions, /(M7/2) <C 1, so that the 
distinction between @q and fi can be disregarded (see 
below). 

Using Eqs. (1) and (3), we arrive at the following for- 
mula for the GLPD dc responsivity: 
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Using Eqs. (2) - (6), we arrive at 
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Assuming that K = 1 and g th = 10 21 cm" 2 s _1 [271 ] . 
at T = 300 K for 0,/2-k = 1-2 THz we obtain 
D* ~ (4.1 - 8.2) x 10 8 cmHz^VW. Setting K = 50, 
we arrive at D* ~ (1.7 - 3.4) x 10 9 cm-Hz 1 / 2 ^. Due 
to a significant decrease in the thermogeneration rate at 
lower temperatures, the detectivity markedly increases 
with decreasing temperature. Indeed, at T = 77 K, set- 
ting g th = 10 13 cm- 2 s- 1 [13], we obtain D* ~ (1.7 - 
3.4) x 10 13 cm.Hz^VW. 

Figure 2 shows the dependences of the dc responsivity 
and detectivity on the number of GLs K calculated for 
n/2ir = 2 THz using Eqs. (5) and (7) for GLPDs with 
sufficiently large / in which the tunneling is weak (see 
below) . 

Equation (7) shows that the GLPD detectivity de- 
creases with increasing photon frequency. This is be- 
cause the spectral dependence of the GLPD detectivity 
is determined by that of the responsivity. As follows 
from Eqs. (6) and (7), the detectivity drops at elevated 
bias voltages when the interband tunneling prevails over 
their thermogeneration of the electron and hole pairs. 
According to Eqs. (4) and (7), the GLPD responsivity 
is independent of the length of the i-section 21 and the 
bias voltage V (at least if the latter is not too small), 
whereas the detectivity increases with increasing / and 
decreases with V. This is because the component of the 
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FIG. 2: Responsivity Ro and detectivity D* versus number 
of GLs K for Q./2-k = 2 and 3 THz. 



FIG. 3: Normalized detectivity versus number bias voltage 
for GLDPs with different length of i-section 21. 



dark current associated with the thermogeneration rate 
increases linearly with increasing I. At the same time, 
the tunneling component is proportional to V 3 ^ 2 /VI. As 
a result, from Eq. (7) one can arrive at 



D* oc 



v/i + uy/W 2 



(8) 



where b oc 1/gth- Thus, at elevated bias voltages 
when the tunneling generation surpassed the thermal 
generation, D* oc V~ 3 / 4 . Figure 3 shows the voltage- 
dependence of the detectivity D* tunn with the interband 
tunneling normalized by the detectivity D* without tun- 
neling calculated for different values of the i-section 
length 21. 

Since the value of D* is determined by both I and V, 
the latter values should be properly chosen. 

There is also a limitation associated with the neces- 
sity to satisfy the condition 21 < Ir, where Ir is the 
recombination length. Otherwise, the recombination in 
the i-section can become essential resulting in degrading 
of the GLPD performance. The recombination length 
can be defined as Ir — Ly, where Ly —< v > tr, tr is 
the recombination time, and < v > is the average drift 
velocity in the i-region. Assuming that for the electron 
and hole densities close to that in the intrinsic graphene 
at T = 300 K one can put tjj = 5x 10" 10 s. Setting for 
sufficiently large bias < v >= 5 x 10 7 cm/s [Hj], one can 
find Ly ~ 250 /an. 



III. FEATURES OF GLPD WITH 
ELECTRICALLY-INDUCED JUNCTION 

The main distinctions of the detectors with doped and 
electrically induced p-i-n-junction are that the Fermi en- 
ergies of electrons and holes in the latter depend on the 
gate voltages V p < in the p-section and V n > in 
the n- section [see Fig. 1(b)] and that these energies and, 
hence, the heights of the barriers confining electrons and 
holes in the pertinent section are different in different 



GLs. This is due to the screening of the transverse elec- 
tric field created by the gate voltages by the GL charges. 
This screening results in a marked decrease in the barrier 
heights of the GLs located in the MGL-structure depth 
(with large indices k) in comparison with the GLs near 
the top. In the following, we set V n = —V v = V g . The 
barriers in question effectively prevent the electron and 
hole injection from the contacts under the reverse bias 
if the barrier heights are sufficiently large. However, the 
electron and hole injection into the GLs with large k leads 
to a significant increase in the dark current (in addition 
to the current of the electrons and hole generated in the 
i-section). This can substantially deteriorate the GLPD 
detectivity. As for the responsivity, it can still be calcu- 
lated using Eqs. (4) or (5). To preserve sufficiently high 
detectivity, the number of GLs should not be too large. 
Indeed, taking into account the contributions of the cur- 
rent injected from the p- and n-regions to the net dark 
current, the detectivity of a GLPD with the electrically- 
induced p-i-n junction can be presented as 
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where oc exp(—s { p'/kBT) is the electron current in- 
jected from the p-region and the hole current injected 

(k) 

from the n-region, e F is the Fermi energy of holes 
(electrons) in the p-region (n-region) of the fc-th GL 
(l<k<K). 

The injected currents are small in comparison with the 
current of the electrons and holes thermogenerated in the 
i-region if 
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Here So is the electron and hole density in the intrinsic 
graphene (at given temperature). 



Considering Eqs. (A5) and (A6) from the Appendix, 
the latter imposes the following limitation on the num- 
ber of GLs in GLDs with the electrically-induced p-i-n 
junctions and the gate voltage V g : 



K < K r 



oc V 

k B Tln(v w '£>o/2irlg th ) 
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The quantity 7 is defined in the Appendix. For 7 = 
0.07 — 0.12, assuming that the Fermi energy in the top- 
most GL e<P = e T F = 100 mcV, T = 300 K, E - 
8 x 10 10 cm" 2 , I = 10 jum, and g th ~ 10 21 cm^s -1 , 
we obtain K max ~ 34 — 58, i.e., fairly large. A decrease 
in I results in smaller K max due to a substantial increase 
in the tunneling current. For instance, if I = 0.1 /Ltm fTlj . 
we obtain K max ~ 7 - 10. 

Thus, by applying gate voltages one can form the p- 
and n-sections with sufficiently high densities and large 
Fermi energies of holes and electrons in GL-structures 
with a large number of GLs. 



IV. DYNAMICAL RESPONSE 

To describe the dynamic response to the signals modu- 
lated with the frequency ui <C fi, one needs to find the ac 
components of the net density of photogenerated electron 
and holes, <5£~ and <5£+, respectively. These components 
are governed by the following equation: 
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As in above, 5 I L 
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In the case of ballistic transport of electrons and holes 
across the i-section, < v >~ vw- If the electron and 
hole transport is substantially affected by quasi-elastic 
scattering, so that the momentum distributions of elec- 
trons and holes are virtually semi-isotropic, one can put 
< v >~ v w /2 (see also, Ref. [28j|). Solving Eq. (12) with 
the boundary conditions 5T^\ x =±i = 0, we obtain 




exp[iuj(l ± x)f < v >] — 1 



SI,, 



(13) 



[13) and considering the Ramo-Shockley the- 
orem 1281. 12p applied to the case of specific contact ge- 
ometry [30j the ac current induced in the side contacts 
(terminal current), can be presented as 
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FIG. 4: Dependence of responsivity (normalized) on modula- 
tion frequency u)/2tt for different Tt. 



Here we have introduced the characteristic transit time 
r t = 1 1 < v >. The feature of the contact geometry (the 
blade-like contacts) was accounted for by using the form- 
factor <?(£) = 1/ny/l — £ 2 [3l|. This is valid because the 
thickness of the side contacts and multiple-GL structure 
under consideration Kd <C I even at rather large numbers 
of GLs K, where d is the spacing between GLs. Similar 
approach was used previously to analyze the dynamic 
response of the lateral p — n junction photodiodes made 
of the standard semiconductors [12] and the graphene 
tunneling transit-time THz oscillator 22]. Integrating in 
Eq. (14), we arrive at the following: 
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where J7n(£) is the Bessel function. Equation (15) yields 
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Using Eq. (16), for the frequency dependent responsivity 
= \SJP hoto \/2lSI LJ , we obtain 



Ru = K 
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Figure 3 shows the responsivity R u normalized by 
its dc value versus the modulation frequency ui calcu- 
lated using Eq. (17) for different values of the transit 
time Tt oc I. At < v >= vy//2 ~ 5 x 10 7 cm/s, the 
range r t = 10 — 30 ps corresponds to the length of 
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the i-section 2Z = 10 — 30 /im. As seen from Fig. 3, 
Rlj/Rq = 1/V2 when the modulation cut-off frequency 
ft = uj t /2n ~ 8.3 - 24.9 GHz. Naturally, in GLPDs with 
shorter i-sections, ft can be markedly larger (although, 
at the expense of a substantial increase in the tunneling 
current). Due to this, GLPDs can be used as ultrafast 
THz and 1R photodetectors. 



V. COMPARISON WITH QWIPS AND QDIPS 

Now we compare the responsivity and detectivity of 
GLPDs with K GLs calculated above and those of 
QWIPs (properly coupled with the incident THz or IR 
radiation) and QDlPs with the same number of QWs. 
The fraction of the absorbed photon flux in one QW 
I3 { ^ W) ~ u^ W) T,f W \ where a^ w ) is the cross-section 
of the photon absorption due to the intersubband transi- 
tions and the donor sheet density. Setting the usual 



values a. 



(QW) 
a. 



2 x 10- 15 cm 2 and H^ W) 



one obtains /3q* W ^ — 0.002. This value is one order 
of magnitude smaller that 0q. Hence, one can neglect 
the attenuation of radiation in QWIPs with K < 100 
(whereas in GLPDs it can be essential) . In such a case, 
the responsivity of QWIPs is independent of the number 
of QWs K in the QWIP structure [H and given by 



K 



Ml p c 



(18) 



where p c is the so-called capture probability which relates 
to the QWIP gain G^ QW ~> as G^ QW ^> = (l-p c )/Kp c ~ 
1/Kp c (33|. Using Eqs. (4) and (18), we obtain 
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39), 



Assuming that K = 50 - 100, (K* ~ 30 - 
0.002 and p c = 0.1, we obtain R /R { QW) -70-90 
The ratio of the detectivities can be presented as 



D* 



Wn^P~c K* 

D*(QW) ~ R (QW) K 
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(20) 



As can be extracted from Ref. [26|], g t h oc 
cxp(—hui /kBT) , where hui ~ 0.02 eV is 
the energy of optical phonon in GLs, whereas 
g^ W) cx eM-e {QW) /k B T), where e« w ) is the 
QW ionization energy (Ml > e {Qw ^), Eq. (17) yields 



D* 



D*(QW) R (QW) 
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H(u>o — fi) 
2k B T 
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Different dependences of the responsivities and detectivi- 
ties of GLDs and QWIPs on K is due different directions 



of the dark current and photocurrent: parallel to the GL 
plane in the former case and perpendicular to the QW 
plane in the latter case. The product of the factors in 
the right-hand side of Eq. (21) except the last one can 
be on the order of unity. This is because a large ra- 
tio 2/?n/ Pq^ W ^ can be compensated by relatively small 
capture parameter p c . However, the exponential factor 
in Eq. (21) is large in the THz range: £1 < u>o, i.e., 
at tt/2ir < 50 THz. In particular, at T = 300 K and 
0/27T = 5 — 10 THz, the exponential factor in question 
is about 24 - 36. 

The GLPD responsivity and detectivity can also 
markedly exceed those of QDIPs (for which formulas sim- 
ilar to Eqs. (21) and (22) can be used) despite lower 
capture probability and thermoexcitation rate in QDIPs 
in comparison with QWIPs (see, for instance, Refs. [H, 
IH, [36|. The ratios D*/D*^ w ^ and D*/D*^ D ^ dra- 
matically increases with decreasing Ml and T. This 
is attributed to the fact that the spectral dependence 
of the GLPDs detectivity is similar to the spectral de- 
pendence of the responsivity (D* cx /3q/MY), whereas 
the GLPDs intended for the photodetection in different 
spectral ranges exhibit the same dark current (due to 
the gapless energy spectrum). In contrast, in QWIPs, 
QDIPs, and some other photodetectors (see, for exam- 
ple, Refs. d,[I3|) the transition to lower photon frequency 
requires to utilize the structures with lower ionization en- 
ergy and, hence, exponentially higher dark current. The 
latter leads to quite different spectral dependencies of 
the GLPD detectivity and the detectivity of QWIPs and 
QDIPs (which drops with increasing photon freqiency). 

GLPDs can surpuss QWIPs and QDIPs in responsivity 
even at relatively high photon frequencies. For example, 
for a GLPD with K = 25 at 0/2tt = 75 THz (wave- 
length A about 4 /im), we arrive at Ro ~ 3 A/W. This 

value is three times larger than R^ D ^ obtained experi- 
mentally for a QDIP with 25 InAs QD layers [13] ■ As for 
the detectivity, comparing a GLPD with K = 70 oper- 
ating at T = 300 K and Q/2n ~ 15 THz (A ~ 20 (im) 
and a QDIP with 70 QD layers [H], we obtain D* ~ 
2 x 10 s cm-Hz^/W and D*^ D ^> ~ 10 7 cm-Hz^/W, 
respectively. 

GLPDs can surpass also photodetectors on narrow-gap 
and gapless bulk semiconductors like HgCdTe (BSPDs). 
Apart from advantages associated with potentially sim- 
pler fabrication, GLPDs might exhibit higher detectivity 
(compare the data above and those from Ref. 39]). This 
can be attributed to relatively low thermogeneration rate 
in GLPDs compared to BSPDs with very narrow or zc- 
roth gap. The point is that thermogeneration rate in 
GLPDs at room temperatures is primarily due absorp- 
tion of optical phonons [27j which have fairly large en- 
ergy huiQ (the Auger processes in GLs are forbidden [40j), 
whereas this rate in BSPDs is essentially determined by 
the Auger processes which are strong |39j . 
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VI. DISCUSSION AND CONCLUSIONS 

Analyzing Eqs. (4) - (7), we accepted that /(Ml) -C 1 
and, therefore, disregarded the frequency dependence of 
the absorption coefficient /3q associated with the popu- 
lation of the low energy states by electrons and holes. 
For a rough estimate, the value of distribution function 
/(0) can be presented as /(0) < /o(0)£/£ , where S is 
the electron and hole density in the i-section under the 
reverse bias. This density, in turn, can be estimated as 
£ = 2g t hT t = gth^l/ < v >. Using the same parame- 
ters as above, for I = 10 /xm and T = 30 K, we obtain 
£/£ = 0.25, so that /(0) < 0.125 and 0n|n^o > 0.75/3. 
It implies that the numerical data obtained above for the 
frequencies at the lower end of the THz range might be 
slightly overestimated. However, since the electron-hole 
system in the i-section can be pronouncedly heated by the 
electric field, the actual values of f(Ml/2) can be smaller 
than in the latter estimate and, hence, /3q can be rather 
close to 0. The electron and hole heating in intrinsic 
GLs under the electric field was studied recently l28l . l4~fl | . 
In contrast to the cases considered in Refs. [28|, [42|, the 



finiteness of the transit time of electrons and holes in 
the i-section of the GL-structures under consideration 
strongly affects the electron and hole heating. There- 
fore, more careful calculation of (3q at relatively low f2 
requires separate consideration. 

In summary, we proposed and evaluated GLPDs 
multiple-GL p-i-n structures. It was shown that GLPDs 
can exhibit high responsivity and detectivity in the THz 
and IR ranges at room temperatures. Due to relatively 
high quantum efficiency and low thermogeneration rate, 
the GLPD responsivity and detectivity can substantially 
exceed those of other photodetectors. 
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Appendix. Effect of vertical screening in multiple 
GL-structures 



As was pointed out above, the thickness of multiple 
GL-structures even with rather large number of GLs K 
is in reality small in comparison with the lateral sizes 
of the device, namely, the lengths of all the section and, 
hence, the gates. Owing to this, the distribution of the dc 
electric potential tp = 2ip /V g normalized by V g /2 (where 
V g = V„ = —V p ) in the direction perpendicular to the GL 
plane (corresponding to the axis z) can be found from the 
one-dimensional Poisson equation: 



8?re 



K 



eeV t 



k=l 



Here Eq^ is the electron (hole) density in the fc-th GL in 
the n-section (p-section), d is the spacing between GLs, 
and 8{z) is the Dirac delta-function. Considering that 

= (4) 2 MH = (e 2 ^ 2 /4^ 2 4)^ 2 | z=fed , where 
e k F is the Fermi energy in the n-section (p-section) of 



the A;-th GL, and replacing the summation in Eq. (Al) 
by integration (that is valid if K is not too small), we 
reduce Eq. (Al) to the following: 



d 2 ^ 
Ih 2 



(A2) 



with the characteristic screening length L s = 
hvp \/as d/2e 3 V g oc V^ 1 ^ 2 . One can assume that 
il>\z=o = 1 + 2Wg(d-i(j/dz)\ z=0 and ^\ z=oc = (as well as 
(dtp/dz)\ z=oc = 0), where W g is the thickness of the layer 
separating the multiple GL-structure and the gates. 
Solving Eq. (A2) with the latter boundary conditions, 
we arrive at 



1 



(C + z/V&L s f ' 
where C satisfies the following equation: 



(A3) 



C 3 -C 2 = (4Wg/V6L s ), (A4) 
Since in reality AW g ^> \/6L s , from Eq. (A4) one obtains 

C ~ (AW g /V6L s )^ 3 oc Vg 1/6 Setting d = 0.35 nm, W g = 
10 nm, ae = 4, and V g — 2 V, one can obtain L 8 ~ 0.3 nm, 
(4Wg/VQL s ) ~ 53 and C ~ 3.75. 
Equation (A3) yields 



t P — 



eV a 



2[C+(k-l)d/V6L s ] 2 



(A5) 



Here = e^P = eV g /2C oc V g 5 ^ e is the Fermi energy 
of electrons in the topmost GL in the n-section (holes in 
the p-section) and 



a« = [l + (fc-l) 7 ]- 2 , 



(A6) 



where 7 = d/V6L s C oc V^ 3 At the above parameters 
(in particular, W g = 10 — 50 nm and V g = 2V) 7 = 
d/V6L s C ~ 0.07-0.12. 



